MD simulations indicate Omicron P132H of SARS-CoV-2 Mpro is a potential allosteric mutant involved in modulating the dynamics of catalytic site entry loop.

The SARS-CoV-2 main protease Mpro, essential for viral replication is an important drug target. It plays a critical role in processing viral polyproteins necessary for viral replication assembly. One of the predominant SARS-CoV-2 Mpro mutations of Omicron variant is Pro132His. Structurally, this mutation site is located ∼22 Å away from the catalytic site. The solved crystal structure of this mutant in complex with inhibitors as well as its reported catalytic efficiency did not show any difference with respect to the wild type. Thus, the mutation was concluded to be non-allosteric. Based on microsecond long MD simulation of the Pro132His mutant and wild type, we show that Pro132His mutation affects the conformational equilibrium with more population of conformational substates having open catalytic site, modulated by the dynamics of the catalytic site entry loop, implying the allosteric nature of this mutation. The structural analysis indicates that rearrangement of hydrogen bonds between His132 and adjacent residues enhances the dynamics of the linker, which in turn is augmented by the inherent dynamic flexibility of the catalytic pocket entry site due to the presence of charged residues. The altered dynamics leading to loss of secondary structures corroborate well with the reported compromised thermal stability.

Physiological models to study the effect of molecular crowding on multi-drug bound proteins: insights from SARS-CoV-2 main protease.

Molecular Dynamics simulations are often used in drug design. However, such simulations do not account for the physiological environment of the receptor; hence overlook its impact on biomolecular interactions. To address this lacuna, we identified three objectives to pursue - develop models of physiological environment, study a drug-receptor complex in such environments, and identify methods to analyze these complicated simulations. Two novel physiological models were developed and studied. The first, called 'm10', comprises of 10 of the most abundant cytoplasmic metabolites at physiological concentrations. The second, called 'phy', supplements m10 with an additional crowder protein to elicit macromolecular crowding effect. The main protease (Mpro) of SARS-CoV-2, being essential for viral replication, is an attractive drug target for COVID-19. Hence, we chose Mpro docked with multiple drugs as our model drug-receptor system. With a plethora of compounds, physiological systems can be exceedingly large and complex. A novel Spark-based software (SparkTraj) was developed to rapidly analyze non-specific contacts and water interactions. Our study shows that crowding enhances the difference in the dynamics of apo- vs drug-bound complexes. Metabolites, at times as a cluster, were seen interacting with the protease, drugs, and binding sites in drug-free receptor. Except one that crawled to an adjacent pocket in phy, the drugs remained in their respective pockets in all simulations. Given these observations, we hope that the models and approach presented here would help the optimization, evaluation, and selection of potential drugs. Generic biomolecular dynamics could also benefit from such models and tools.Communicated by Ramaswamy H. Sarma.

Targeting allosteric pockets of SARS-CoV-2 main protease Mpro.

Repurposing of antivirals is an attractive therapeutic option for the treatment of COVID-19. Main protease (Mpro), also called 3 C-like protease (3CLpro) is a key protease of SARS-CoV-2 involved in viral replication, and is a promising drug target for antivirals. A major challenge to test the efficacy of antivirals is the conformational plasticity of Mpro and its future mutation prone flexibility. Suitable choice of drugs in catalytic and allosteric pockets appear to be essential for combination therapy. Current study, based on docking and extensive set of MD simulations, finds the combination of Elbasvir, Glecaprevir and Ritonavir to be a viable candidate for further experimental drug testing/pharmacophore design for Mpro.Communicated by Ramaswamy H. Sarma.

The p53-Mdm2 interaction and the E3 ligase activity of Mdm2/Mdm4 are conserved from lampreys to humans.

The extant jawless vertebrates, represented by lampreys and hagfish, are the oldest group of vertebrates and provide an interesting genomic evolutionary pivot point between invertebrates and jawed vertebrates. Through genome analysis of one of these jawless vertebrates, the Japanese lamprey (Lethenteron japonicum), we identified all three members of the important p53 transcription factor family--Tp53, Tp63, and Tp73--as well as the Mdm2 and Mdm4 genes. These genes and their products are significant cellular regulators in human cancer, and further examination of their roles in this most distant vertebrate relative sheds light on their origin and coevolution. Their important role in response to DNA damage has been highlighted by the discovery of multiple copies of the Tp53 gene in elephants. Expression of lamprey p53, Mdm2, and Mdm4 proteins in mammalian cells reveals that the p53-Mdm2 interaction and the Mdm2/Mdm4 E3 ligase activity existed in the common ancestor of vertebrates and have been conserved for >500 million years of vertebrate evolution. Lamprey Mdm2 degrades human p53 with great efficiency, but this interaction is not blocked by currently available small molecule inhibitors of the human HDM2 protein, suggesting utility of lamprey Mdm2 in the study of the human p53 signaling pathway.

Structural basis for microtubule recognition by the human kinetochore Ska complex.

The ability of kinetochores (KTs) to maintain stable attachments to dynamic microtubule structures ('straight' during microtubule polymerization and 'curved' during microtubule depolymerization) is an essential requirement for accurate chromosome segregation. Here we show that the kinetochore-associated Ska complex interacts with tubulin monomers via the carboxy-terminal winged-helix domain of Ska1, providing the structural basis for the ability to bind both straight and curved microtubule structures. This contrasts with the Ndc80 complex, which binds straight microtubules by recognizing the dimeric interface of tubulin. The Ska1 microtubule-binding domain interacts with tubulins using multiple contact sites that allow the Ska complex to bind microtubules in multiple modes. Disrupting either the flexibility or the tubulin contact sites of the Ska1 microtubule-binding domain perturbs normal mitotic progression, explaining the critical role of the Ska complex in maintaining a firm grip on dynamic microtubules.

Conservation of all three p53 family members and Mdm2 and Mdm4 in the cartilaginous fish.

Analysis of the genome of the elephant shark (Callorhinchus milii), a member of the cartilaginous fishes (Class Chondrichthyes), reveals that it encodes all three members of the p53 gene family, p53, p63 and p73, each with clear homology to the equivalent gene in bony vertebrates (Class Osteichthyes). Thus, the gene duplication events that lead to the presence of three family members in the vertebrates dates to before the Silurian era. It also encodes Mdm2 and Mdm4 genes but does not encode the p19(Arf) gene. Detailed comparison of the amino acid sequences of these proteins in the vertebrates reveals that they are evolving at highly distinctive rates, and this variation occurs not only between the three family members but extends to distinct domains in each protein.

A novel p53 phosphorylation site within the MDM2 ubiquitination signal: II. a model in which phosphorylation at SER269 induces a mutant conformation to p53.

The p53 DNA-binding domain harbors a conformationally flexible multiprotein binding site that regulates p53 ubiquitination. A novel phosphorylation site exists within this region at Ser(269), whose phosphomimetic mutation inactivates p53. The phosphomimetic p53 (S269D) exhibits characteristics of mutant p53: stable binding to Hsp70 in vivo, elevated ubiquitination in vivo, inactivity in DNA binding and transcription, increased thermoinstability using thermal shift assays, and λ(max) of intrinsic tryptophan fluorescence at 403 nm rather than 346 nm, characteristic of wild type p53. These data indicate that p53 conformational stability is regulated by a phosphoacceptor site within an exposed flexible surface loop and that this can be destabilized by phosphorylation. To test whether other motifs within p53 have similarly evolved, we analyzed the effect of Ser(215) mutation on p53 function because Ser(215) is another inactivating phosphorylation site in the conformationally flexible PAb240 epitope. The p53(S215D) protein is inactive like p53(S269D), whereas p53(S215A) is as active as p53(S269A). However, the double mutant p53(S215A/S269A) was transcriptionally inactive and more thermally unstable than either individual Ser-Ala loop mutant. Molecular dynamics simulations suggest that (i) solvation of phospho-Ser(215) and phospho-Ser(269) by positive charged residues or solvent water leads to local unfolding, which is accompanied by local destabilization of the N-terminal loop and global destabilization of p53, and (ii) the double alanine 215/269 mutation disrupts hydrogen bonding normally stabilized by both Ser(215) and Ser(269). These data indicate that p53 has evolved two serine phosphoacceptor residues within conformationally flexible epitopes that normally stabilize the p53 DNA-binding domain but whose phosphorylation induces a mutant conformation to wild type p53.

Differential binding of p53 and nutlin to MDM2 and MDMX: computational studies.

Half of human tumours have mutated p53 while in the other half, defective signalling pathways block its function. One such defect is the overexpression of the MDM2 and MDMX proteins. This has led to an intense effort to develop inhibitors of p53-MDM2/MDMX interactions. Nutlin is the first such compound described to block p53-MDM2 interactions. Molecular dynamics simulations have been used to explore the differences in binding of p53 and nutlin to MDM2/MDMX. Simulations reveal that p53 has a higher affinity for MDM2 than MDMX, driven by stronger electrostatic interactions. p53 is displaced from MDM2 by nutlin because it is more flexible, thus paying a larger entropic penalty upon sequestration by MDM2. The inherent plasticity of MDM2 is higher than that of MDMX, enabling it to bind both p53 and nutlin. The less flexible MDMX interacts with the more mobile p53 because the peptide can adapt conformationally to dock into MDMX, albeit with a reduced affinity; nutlin, however is rigid and hence can only interact with MDMX with low affinity. Evolutionarily, the higher affinity of MDM2 for p53 may enable MDM2 to bind p53 for longer periods as it shuttles it out of the nucleus; in contrast, MDMX only needs to mask the p53 TA domain. This study enables us to hypothesize gain of function mutations or those that have decreased affinity for nutlin. These conclusions provide insight into future drug design for dual inhibitors of MDM2 and MDMX, both of which are oncoproteins found overexpressed in many cancers.

Differences in the transactivation domains of p53 family members: a computational study.

The N terminal transactivation domain of p53 is regulated by ligases and coactivator proteins. The functional conformation of this region appears to be an alpha helix which is necessary for its appropriate interactions with several proteins including MDM2 and p300. Folding simulation studies have been carried out to examine the propensity and stability of this region and are used to understand the differences between the family members with the ease of helix formation following the order p53 > p73 > p63. It is clear that hydrophobic clusters control the kinetics of helix formation, while electrostatic interactions control the thermodynamic stability of the helix. Differences in these interactions between the family members may partially account for the differential binding to, and regulation by, MDM2 (and MDMX). Phosphorylations of the peptides further modulate the stability of the helix and control associations with partner proteins.

The Mdm2 and p53 genes are conserved in the Arachnids.

The p53 protein and its negative regulator the ubiquitin E3 ligase Mdm2 have been shown to be conserved from the T. adhaerens to man. In common with D. melanogaster and C. elegans, there is a single copy of the p53 gene in T. adhaerens, while in the vertebrates three p53-like genes can be found: p53, p63 and p73. The Mdm2 gene is not present within the fully sequenced and highly annotated genomes of C. elegans and D. melanogaster. However, it is present in Placazoanand the presence of multiple distinct p53 genes in the Sea anemone N. vectensis led us to examine the genomes of other phyla for p53 and Mdm2-like genes. We report here the discovery of an Mdm2-like gene and two distinct p53-like genes in the Arachnid Ioxodes scapularis (Northern Deer Tick). The two predicted Deer Tick p53 proteins are much more highly related to the human p53 protein in sequence than are the fruit fly and nematode proteins. One of the Deer Tick genes encodes a p53 protein that is initiated within the DNA binding domain of p53 and shows remarkable homology to the newly described N-terminally truncated delta isoforms of human and zebrafish p53.

Mdm2 and p53 are highly conserved from placozoans to man.

The p53 protein is the most commonly mutated tumor suppressor gene in man. Understanding of its evolutionary origins have been enhanced by the recent discovery of p53 family genes in the Sea Anemone Nematostella vectensis. This amino acid sequence conservation has been reflected in biological activity since the early p53 proteins, like their human counterparts, are responsible for DNA damage-induced cellular apoptosis, albeit restricted to the germ cell compartment in model organisms such as the nematode and fruit fly. In vertebrates from zebrafish to man the function of p53 is tightly and absolutely constrained by a negative regulator Mdm2. However the Mdm2 gene has not been detected in the genome of the model nematode (C. elegans) and insect (D. melanogaster) species. We have found that the p53 gene and the Mdm2 gene are present in Placozoans, one of the simplest of all free living multi-cellular organisms, implying that both proteins arose much earlier in evolution than previously thought. Detailed sequence analysis shows the exceptional retention of key features of both proteins from man to Placazoan implying that the p53-Mdm2 interaction and its regulation have been conserved from a basal eumetazoan since the pre-cambrian era over 1 billion years ago.

Design of a novel MDM2 binding peptide based on the p53 family.

p53 is a major tumor suppressor protein, that binds to, and is negatively regulated by MDM2. In tumors overexpressing MDM2, p53 function can be rescued through the disruption of the MDM2-p53 interactions by small molecules and peptides. It is known that MDM2 also binds p73 but not p63, the two homologues of p53. We dissect the structural and energetic reasons underlying this discrimination and have identified a peptide that is intrinsically less helical than p53 and yet has a higher affinity for MDM2. The increased disorder has been introduced by localizing a cationic residue in between two anionic residues, imparting a degree of frustration to the system. In addition, the introduction of a bulkier hydrophobic group towards the centre of the peptide enables the peptide to adapt a bound conformation that on the one hand is most strained, and yet enables the peptide to straddle the largest surface of MDM2, amongst all the peptides. Computations also reveal that this peptide is a dual inhibitor, binding also to MDMX. The computed affinity of the new peptide has been validated against MDM2 using fluorescence-based thermal shift assays.

Dimerization of the core domain of the p53 family: a computational study.

Computational models reveal the structural origins of cooperativity in the association of the DNA binding domains (DBD) of p53 (and its two homologues p63 and p73) with consensus DNA. In agreement with experiments they show that cooperativity, as defined by sequential binding of monomers to DNA is strong for p53 and weak for homologues p63 and p73. Computations also suggest that cooperativity can arise from the dimerization of the DBD prior to binding the DNA for all 3 family members. Dimerization between the DBDs is driven by packing interactions originating in residues of helix H1 and loop L3, while DNA binding itself is dominated by local and global electrostatics. Calculations further suggest that low affinity oligomerization of the p53 DBD can precede the oligomerization of the tetramerization domain (TD). During synthesis of multiple chains on the polysome, this may increase fidelity by reducing the possibility of the highly hydrophobic TD from nonspecific aggregation. Mutations have been suggested to test these findings.

Stability of the core domain of p53: insights from computer simulations.

BACKGROUND: The tumour suppressor protein p53 protein has a core domain that binds DNA and is the site for most oncogenic mutations. This domain is quite unstable compared to its homologs p63 and p73. Two key residues in the core domain of p53 (Tyr236, Thr253), have been mutated in-silico, to their equivalent residues in p63 (Phe238 and Ile255) and p73 (Phe238 and Ile255), with subsequent increase in stability of p53. Computational studies have been performed to examine the basis of instability in p53. RESULTS: Molecular dynamics simulations suggest that mutations in p53 lead to increased conformational sampling of the phase space which stabilizes the system entropically. In contrast, reverse mutations, where p63 and p73 were mutated by replacing the Phe238 and Ile255 by Tyr and Thr respectively (as in p53), showed reduced conformational sampling although the change for p63 was much smaller than that for p73. Barriers to the rotation of sidechains containing aromatic rings at the core of the proteins were reduced several-fold when p53 was mutated; in contrast they increased when p73 was mutated and decreased by a small amount in p63. The rate of ring flipping of a Tyrosine residue at the boundary of two domains can be correlated with the change in stability, with implications for possible pathways of entry of agents that induce unfolding. CONCLUSION: A double mutation at the core of the DNA binding domain of p53 leads to enhanced stability by increasing the softness of the protein. A change from a highly directional polar interaction of the core residues Tyr236 and Thr253 to a non-directional apolar interaction between Phe and Ile respectively may enable the system to adapt more easily and thus increase its robustness to structural perturbations, giving it increased stability. This leads to enhanced conformational sampling which in turn is associated with an increased "softness" of the protein core. However the system seems to become more rigid at the periphery. The success of this methodology in reproducing the experimental trends in the stability of p53 suggests that it has the potential to complement structural studies for rapidly estimating changes in stability upon mutations and could be an additional tool in the design of specific classes of proteins.